Fermion pair production in QED and the backreaction problem 
in (l+l)-dimensional boost-invariant coordinates revisited 



Bogdan Mihaila,^'0 John F. Dawson,^'0 and Fred Cooper'^' "'•[|| 

Los Alamos National Laboratory, Los Alamos, NM 87545 
^ Department of Physics, University of New Hampshire, Durham, NH 03824 
'^National Science Foundation, 4201 Wilson Blvd., Arlington, VA 22230 
^ Santa Fe Institute, Santa Fe, NM 87501 

We study two difFerent initial conditions for fermions for the problem of pair production of fermions 
coupled to a classical electromagnetic field with backreaction in (1+1) boost-invariant coordinates. 
Both of these conditions are consistent with fermions initially in a vacuum state. We present results 
for the proper time evolution of the electric field E, the current J, the matter energy density e, 
and the pressure p as a function of the proper time for these two cases. We also determine the 
interpolating number density as a function of the proper time. We find that when we use a "first 
order adiabatic" vacuum initial condition or a "free field" initial condition for the fermion field, 
we obtain essentially similar behavior for physically measurable quantities. The second method is 
computationally simpler, it is twice as fast and involves half the storage required by the first method. 
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I. INTRODUCTION 

Particle production from strong fields has a long his- 
tory starting with Schwinger's classic paper A de- 
tailed history of this subject can be found in two re- 
cent reviews 0, Q. One of the many applications of 
pair production has been as a model for particle produc- 
tion in the central rapidity region following a rclativis- 
tic heavy ion collision. Following such a collision, there 
is experimental evidence that the production of parti- 
cles is "boost-invariant" [1, Q which leads to measurable 
quantities such as energy densities being functions of the 
fluid proper time alone. The initial conditions we want 
to study for this problem are that the number of pairs 
starts out zero and that the initial induced current in 
the Maxwell (backreaction) equation for the electric field 
is also zero. We then want to study the proper time 
evolution of the expectation value of the energy density, 
pressure, and current of the produced particles and the 
evolution of the electric field. 

In an earlier paper on this topic one particular set 
of initial conditions consistent with having no pairs of 
particles produced before the collision at initial proper 
time, T = To, led to the need for doubling the number 
of fermion solutions in order to start with zero induced 
current in the Maxwell equation for the electric field. In 
the paper by Cooper et al. @ , two sets of solutions for 
the second-order squared Dirac equation were used in or- 
der to satisfy the desired initial conditions of having zero 
initial current. Similar results were presented in a 1992 
paper by Kluger et al. for the Cartesian case. In those 
papers, the initial conditions were taken to correspond to 
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a first-order adiabatic approximation to the second-order 
Dirac equation, which forced them to average over two 
different solutions of the Dirac equation so that the cur- 
rent vanished at r = tq. In the present work, we consider 
a slightly different initial state, namely approximate free 
fields for the fermions, which automatically leads to a 
zero current at r = tq. This initial condition was used 
earlier by Cooper and Savage [1] in their study of the 
dynamics of the chiral phase transition in the (2-f 1) di- 
mensional Gross-Neveu model. The free field initial con- 
dition does not require doubling the number of solutions 
as did the adiabatic choice. We compare the evolution 
of the problem for both initial conditions and show that 
at short to moderate times they are equivalent and are 
slightly different at very late times. 

This semi-classical approximation to the initial value 
QED problem describes the fermions as a quantum field 
but treats the electric field classically. The current used 
in Maxwell's equation is calculated using the vacuum ex- 
pectation value of the quantum Dirac current. As dis- 
cussed in previous papers Q , this approximation is equiv- 
alent to the first term in a large-N approximation to N- 
QED where there are N flavors of fermions present. 

The method we use for numerically solving this prob- 
lem is a shooting method to numerically step out solu- 
tions of the equations from initial conditions. An adia- 
batic analysis of the form of the solutions is used to de- 
termine the behavior of the solutions at large momentum 
and to isolate divergences and perform rcnormalization 
as well as to choose appropriate initial states for fermions 
that are appropriate vacuum states. 

We study here the problem in (1+1) boost-invariant 
coordinates. This kinematic situation is related to the 
kinematics of the early phase of plasma evolution fol- 
lowing a relativistic heavy ion collision with the electric 
field a simplification for the semiclassical chromoelectric 
field expected to be produced in that situation. We study 
(1+1) dimensions for simplicity here, where charge renor- 
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malization is finite. However, the same methods of solu- 
tion used here can be appHed to the case of (3+1) dimen- 
sions in both Cartesian and boost-invariant coordinates. 
We will present results for (3+1) dimensions for QED 
and QCD elsewhere. 

This paper is organized as follows: In Sec. |TT] we 
review briefly the equations we will need to solve for 
QED in (1+1) dimensional boost-invariant coordinates 
(for a detailed derivation, sec e.g. Rcfs. i and|I3). In 
Sees, mil and IIVI wc discuss the backreaction equation 
and the adiabatic expansion and charge renormalization, 
whereas in Sec.|V]we review the calculation of the energy- 
momentum tensor. The two types of initial conditions are 
introduced in Sec. IVII We present results of our numeri- 
cal simulations in Sec. IVIII and conclude in Sec. IVIIII 



II. NOTATION AND EQUATIONS 

In our simplified kinematics, in (1+1) dimensions, we 
choose the longitudinal axis of the collision to be the z- 
axis. Then, in the Cartesian frame, we want to solve the 
set of equations: 

{7''[*9,-eA,(0]-m}^.(0 = 0, (2.1) 

where ^ is shorthand for the Cartesian pair (t, z), '0(^) is 
a fermi field satisfying the anti-commutation relation: 



{ V'a(^, V-K^', } = ~ z') 



(2.2) 



and Aa{(,) is a classical field satisfying Maxwell's equa- 
tions: 

da F^^i^) J^iO, F^^O = d^'A^iO-d^A'^iO- (2.3) 
The current is given by: 

^''^'^^ - (2.4) 



The 7-matrices satisfy {7^,7''} = 277°''' and are given 
by: 



by: 



Boost-invariant coordinates x'^ — [t^t]) are defined 



t ~ T cosh 77 , z = T sinh rj . 



(2.5) 



The connection between the Cartesian frame (d^"), and 
the boost-invariant frame (dx^) is described by a vierbein 
matrix y^(a;), given by: 

dr = V%{x) dx^ , = V\{x) da , (2.6) 

ya j-^-j _ _ f cosh 77, T sinh 77^ 



dxf^ \^ sinh 77, TCOshr]J ' 



and its inverse: 

dx^ = V'aix) dC , da = V'aix) 9^ , (2.7) 

,s _ dx^ _ / cosh77, — sinh77^ 
° 9^" I— sinh7//r, COSI177/T 

The 7-matrices in this frame are denoted by a tilde: 
^^(x) = y^a(x)7", and satisfy: 

{ rix),rix) } = 2g^^x) , g^'^ix) = diag( 1, -l/r^ ) . 

(2.8) 

Dirac's equation (|2.1[) becomes in this frame: 



[r{x)itdf,-eAf,{x))~m]iix) = 0, (2.9) 

where the fermi field obeys the anti-commutation rela- 
tion: 

{ i'air, 77), 4, (t, rj') } = iv) S{V - v')/r • (2.10) 

However, it is much simpler to make a similarity trans- 
formation to a system of coordinates where the vierbein 
becomes diagonal [l^. In this rotated system, the 7- 
matrices are denoted by a bar: 

S-\v)rix)Si^)=rir), (2.11) 

where the matrix S{r]) is given by: 

S'(?7) =exp[77 7V2] = cosh(77/2) + 7^ sinh(77/2) . (2.12) 

The 7^(t) matrices are given explicitly by: 

7- = 7O , 7"(r) = 7Vr , (2.13) 

So if we define a new fermi field (f>{x) by: 

^{x)^S{ii)${x)/V^, (2.14) 

Dirac's equation (|2.9p becomes: 

[z7^(t)V^ -m]0(T, 77)/ V^=0, (2.15) 



where V, 



df, + H^(a;) + ieA^{x), with Ilf,{x) 



S ^(r/) {dfj,S{ri) ), is the covariant derivative. Here 0(x) 
obeys the simpler anti-commutation relation: 

{ kir, V), $1 (r, v') } = -5(77 - 77') . (2.16) 

For our case, the only non- vanishing H^(.t) is for H^ = 
7^/2. In the boost-invariant frame, we work in the tem- 
poral gauge and choose Afj_{x) ~ (0, — A(t)). That is 
A{t) as the negative of the covariant component in the 
boost-invariant frame. So (|2.15p simplifies to: 

{7709^+73 [td,j + eA{T)]/T-m}${T,r]) = 0. (2.17) 



by: 



We now expand the field (/)(t, 7^) in a fourier series given 



/+00 
[dk] Ai'^ e^>^^ ^i^\r) , (2.18) 
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where we have introduced the notation [dk] = dA:/(27r)]. 
Here, A^^^^ are mode operators and (p^j^"^ (r) arc two inde- 
pendent mode functions satisfying the equation: 

[i7°a.-73 7rfc(r)-m] cf>i^\T) ^ , (2.19) 

where 

^fc(T) = - [fc-eA(T)] . (2.20) 
r 

It is now useful to add and subtract the upper and lower 
components of the spinor i/ij.'^'' by writing: 

^i'\r)^UFi'\r), with^=i=(^J . (2.21) 

Here t/^ = U^^ = U'^ . Then f'i,^\t) satisfies an equa- 
tion of Hamiltonian form: 

ldrFlt\T)=H{T)Fl^\T), (2.22) 

with 

Here Kfe(T) is a vector defined in an abstract space TZ 
with imit vectors (ei, 62, e^) and given by: 

Kfe(T)=mei+7rfe(T)e3. (2.24) 

We can introduce the 2x2 dimensional density matrix 
Pk{T) and a "polarization" vector Pfe(T) in TZ with the 
definitions: 

P^k^{r)^Fi''\T)Fi^^\T) = i(l + Pi')(r).<T). (2.25) 

Then from (|2.22p . the polarization vector P^'^''(r) obeys 
the vector equation of motion: 

a.Pi')(r) = 2Kfe(r)xp(^)(r). (2.26) 

Since H(t) in Eq. (|2.23|) is hermitian, fI,^\t) satisfies a 
conservation equation: 

dAFi^''\r)FP{r)]=Q. (2.27) 

So if we choose the two spinors to be orthonormal at 
T = To, they remain orthonormal for all r. In Sec. IVII 
we show how to do this. So we can assume that these 
spinors are orthonormal and complete for all r: 

Fi^^\T)F[^'\T)^5^,^,, (2.28a) 
Y,Flt\T)Fi^^\r) = l. (2.28b) 

Probability conservation also requires that the polariza- 
tion vector P^'^'' (r) for both of these solutions to remain 



on the unit sphere for all time r. So to summarize, the 
fermi field can be written as: 

V'(r,ry) =5(77)[/F(r,77)/V^, (2.29) 

where the field Fir, rj) obeys the anti-commutation rela- 
tion: 

{ F„(r, ,/), Fl, (r, r/) } = 6^^^, <5(r/ - 77') ■ (2.30) 

and is expanded in terms of the spinors Fj,^'^ (r) which 
satisfy (|2^ : 

/+00 
[dk] Y,Ai'K^'^^Fl:^\r). (2.31) 

A=± 

We can use this orthogonality to invert (|2.3ip to get: 

/+00 
d^e-^''^Ft'^\T)F{r,rj), (2.32) 
-00 

for any time r. Using ()2.30p . we then find that the mode 
operators A^j^l obey the anti-commutation relation: 

{A[^\Ai^'^^} = {27r)S,^,,Sik-k'). (2.33) 

It is traditional to define separate positive and negative 
energy operators by setting: 

A[+^=a,, and 4"' = • (2.34) 

We choose our initial state to be the vacuum with no 
particle or anti-particle present. Then: 

ajO)=0, and 6;^. |0)=0. (2.35) 
This means that: 

( [ 4'^ ^ 4''^ ] > = - (2^) A -5,,,, Sik - k') , (2.36) 
a result we will use in the next section. 



III. MAXWELL'S EQUATION 

Maxwell's equation is given in Cartesian coordinates 
in Eq. (j2.3p with the current given in Eq. (|2.4p . For our 
boost-invariant coordinates, Maxwell's equation reads: 

-^d^[^gF^^'^{x)]^r{x), (3.1) 

where ^/— g = t. Now A^ = (0, — A(r)), so the only 
non-vanishing elements of the field tensor are: 

Fr.^{x) = -Fr,Ax) = -drA{T) = T E{t) (3.2) 

This last equation defines what we call the electric field 
E{t) = {drA{T))/T. Then using the metric g'"'{x) = 
diag( 1, — 1/r^ ), we get: 

F^'"(r) = -F"'^(r) = -E{t)/t , (3.3) 
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and Maxwell's equation becomes: 

drE{T) = - J(t) . (3.4) 
Here we have defined a "reduced" current J(t) by: 



e . . (3-5) 

-:^([0^(^,t),7'0(77,t)]), 

Z T 



Using the field expansion (|2.18p and the expectation 
value (|2.36|) of the mode operators, we find for the re- 
duced current: 



^ ' J ~oo , I , J-oo I 



A=±l • 



A' = ±l 



r' + oo 



/-t-oo 
-OO , 11 



A=±l 



el [dTTk] Pt\^k,r) 



(3.6) 



Here we have used the completeness statement (|2.28bp 
to write the current in terms of positive energy solu- 
tions only. In the last line, we changed integration vari- 
ables from k to TTkir), using dTr^ = dfc/r, and defined 
P(7rfc,r) = Pfc(r). Maxwell's equation p.4p becomes: 

/-(-OO 
[dTTk] Pt\^k,T) . (3.7) 
-OO 

Recall that P3 is the third component of the polarization 
vector in the space TZ. 



where lu = -^/tt^ + and 

1 TTfc^ 



1 TTfc TTfc 

4 w5 



5 TTfc TTfe^ 



(4.3) 



We have suppressed the r dependence here of these quan- 
tities. The dot denotes a partial derivative with respect 
to T. Explicitly, we find: 



Pi = hem 

LU 



P2 = em 



2tj3 



T^k 2 2 

e m 



1 TTfe 

4 u;5 



!_ TTfc TTfc 
4 Cj5 



5 TTkTrl 



(4.4) 



So setting e 



1, Maxwell's equation (|3.7p beeomes: 



+00 



[d^fc] 



1 TTfc 5 TTfc 7T^. 



(4.5) 

All terms odd in ttj. vanish by symmetric integration. 
After integration, Eq. (|4.5p becomes: 



P(r) 



67r m? 



E{t) + J^"^(t) . 



(4.6) 



Here, the first term corresponds to finite charge renor- 
malization in (1-1-1) dimensions and can be brought over 
to the left hand side of the equation. The current J^^^ (t) 
is explicitly finite by power counting and is initially zero. 

An adiabatic expansion of the Dirac equation can also 
be carried out from solutions of the second-order form of 
the Dirac equation. In Section IVl Bl below, we show that 
this gives the same result as in Eqs. (|4.4p . 



ENERGY-MOMENTUM TENSOR 



IV. ADIABATIC EXPANSION 

The large momentum behavior of the solutions of the 
Dirac equation can be obtained by looking at the adia- 
batic expansion of these solutions. Perhaps the simplest 
way to do this is from the polarization equation (|2.26p . 
In order to count powers of time derivatives, we put: 
dr f—^ edr, and set: 



P,{r) = Pf{r) + ePi'\r) + e'P^^\r) 



(4.1) 



Substitution of this into Eq. (|2.26p and equating powers 
of e give the results: 



,(0) 



,(1) 



LU 

Kfc X Kfc 
2a;3 



,(2) _ 3 (Kfc -Kfc) kfc Kfc 



4^5 



(4.2a) 
(4.2b) 

-AAfcK, (4.2c) 



In the boost-invariant coordinate system, the average 
value of the total energy-momentum tensor is given by 
Eqs. (4.1) and (4.2) of Ref.@, and is the sum of two terms: 



T,. = + r^,^" = diag( £,t'V) 



(5.1) 



where 



T-"- = j{[^{x),j^^ix) i^D^^^P{x))]+h.c.) (5.2a) 
1 



1 



r: 



llu'^ — .9^" ^ F^^Faf} + F^a g"^ Fpy 



(5.2b) 



Here = d^-\-ie A^{x) and the subscript notation (/i, v) 
means to symmetrize the term. From our definitions in 
Section [Hi] and Eq. (|3.2p , the field part of the energy- 
momentum tensor is given by: 



Tfi,-'d = diag(i?V2,-r2£;V2). 



(5.3) 



We denote the matter part of the energy-momentum ten- 
sor as: 



T-^"- = diag(e,rV] 



(5.4) 
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For the matter field, we first note that Di, ip{x) = 
S{x)'S/ 1, (f>{x) / where Vi, is the covariant derivative 
defined below Eq. (|2.15p . For the T^r = compo- 
nent, Vo — dr, and using the field expansion (|2.18p . 
Eqs. (EIJ), (11221), and (jOS)) . we find: 

1 r + co 

e{r) = -- / [dk]Y,XTiipi'\r)H{r)] 



— oo 

+ OC 



[dfc]Tr[pW(r)i/(r)] 

> 

/+00 
[dnk]K{TTk)-P^+H^k,r). (5.5) 
-oo 

For the Tr^n = t'^p{t) component, V,, = 9,, — ieA{T) + 
7^/2. Following similar steps to the preceding calcula- 
tion, we find: 

1 r + co 

= - / [dfc]7r,(r)Tr[pi+)(r)a3] 

/ + 00 
[dnk]7TkPt\T^k,T). (5.6) 
-oo 

So from (|5A|) . 

[dTTk] Kink) ■ P(+^ (iTk, r) + — , (5.7a) 

-oo ^ 

/+00 7712 

M ^fcPr'(^fc,r)- — . (5.7b) 
-00 ^ 

The covariant derivative of the energy-momentum ten- 
sor in boost-invariant coordinates is conserved: 

r^^M = + T'^^^T'"' + r;;^r^'" = o . (5.8) 

The Christoffel symbols are defined by: r^^(x) = 
V^a{x) (9^1^°!/ (a;)). In our case, the non-vanishing sym- 
bols are given by: 

r;, = r , r;!^ = = i/r . (5.9) 

So we find that 

dr T^^ T^Vr -f T r'" = , (5.10) 

or 

driT£)+r ^0. (5.11) 



Using the equation of motion (j2.26p and Maxwell's equa- 
tion (|3.7p . one can show that Eq. (jS.lip is automatically 
satisfied. 

Using Eqs. (|4.4p . the adiabatic expansion for the en- 
ergy is given by: 



e2 \ 1 



67rm2 / 2 247rT2 



+ 00 



[dTTk] 2 L 



We recognize the first term as a finite renormalization of 
the charge, the second term as a renormalization of the 
cosmological constant, and the third term as a sum of the 
zero point energies of pairs of particles and anti-particles 
with energy uj{T:k)- We subtract these terms from the 
calculation of the energy and arrive at a finite energy 

£SUb gjypjj |-,y. 



+ 00 



gsu,^ [dTTfc] -K(7rfc).P(7rfc,r)+c^-^ 



For the pressure, the adiabatic expansion gives: 



(5.13) 



E^ 
2 



1 



+00 0_2 
[dTT,]!^ 



67r TO^ / 2 St: 

(5.14) 

Again, the first term renormalizes the charge, the second 
term in canceled by the cosmological constant term and 
the third is the usual pressure. We subtract these terms 
from the pressure to get: 



E^ 



-\-oo 



7^kP3{TTk,T) (5.15) 

w5 8 w7 



Eqs. (|5.13p and (|5.15p are now finite. 



VI. INITIAL CONDITIONS 

The simplest choice of initial conditions is to find ap- 
proximate free-field solutions of Eq. (|2.22p near t = tq. 
This strategy was used in Ref. d, and automatically pro- 
vides a zero current at r = rg. We call this the "one- 
field" method, and is discussed in Section lVI Al below. In 
previous studies of the backreaction problem by Cooper 
et al. adiabatic initial conditions were used which re- 
quired averaging over two different solutions to the Dirac 
equation to obtain an zero current at initial proper time 
Tq. We call this the "two-field" method. We discuss this 
method in Section IVl Bl 



A. One-field method 

At r = To = 1/m, A{tq) = and H{tq) is given by: 

H[r,) = m {^^ . (6.1) 

So at T w Tq, Fk{T) obeys the approximate equation of 
motion: 



Writing 



(5.12) 



«5,Fo;fc(T) = i?(ro)Fo;fc(r) 



i^O;fc(T) =Fo;fce-*"(^-^») 



(6.2) 



(6.3) 
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T 

FIG. 1: (Color online) Proper-time evolution of the electro- 
magnetic fields and current for the one-field and two-field 
methods described in text. Here we choose m = 1, A{tq) — 
and E{tq) = 4. 



We find that uj{tq) = ±wo, where luq ~ m^k'^ + 1. Posi- 
tive frequency solutions given by: 



~(+) _ ^ l ujQ + mk fl\ _ / cos(0fc/2)' 

lsin(0fe/2)^ 



0;fe 



2t^o VC 



(6.4) 



and negative frequency solutions by: 



(-) 



0;fc 



2uJo 



1 



-sin(0,V2)' 
, cos(0fe/2) ^ 



(6.5) 



for these solutions are given by: 



(6.6a) 



cos2(0fc/2) sin(0fc/2) cos(0fc/2) 
^sin(0fc/2) cos(0fe/2) sin\ek/2) 

P^^^^o^W^oy^^) (6-6b) 

sin2(0fc/2) -sin(0fe/2) cos(0fc/2)^ 

^-sin(0fe/2) cos(0fc/2) cos2(0fc/2) 

and are independent of r. The corresponding polariza- 
tion vectors are also independent of r and are given by: 



>(+) 

0:fe 



sin 6'fc ei + cos 9k 63 



Kfe(To) 



UJQ 



p^^^ (6.7) 



The initial spinors are orthogonal and complete: 



(6.8a) 
(6.8b) 



So if we set fI.^\to) — Fq^j^ at r = tq, then the exact 
solutions remain orthogonal and complete for all t and 
(|2.28[) is satisfied. As we have seen in Section lllll only the 
positive energy solutions are needed for the backreaction 
calculation. 

The initial spinors can serve to define a particle num- 
ber operator. Since these initial mode functions form a 
complete set, we can expand the quantum field in terms 
of them: 

/+00 
mY.Al':l{r)e^''^Fl,%{r), (6.9) 
-00 \ 



where ^^^^(r) are mode operators for the Fq^^.i^{t) func- 
tions, which now depend on time. Inverting (j6.9p . we 
find: 



/+00 
dxY,e-^'^^F^,'l;ir)FM,v), (6. 
-00 



10) 



from which we obtain the equal time anti-commutation 
relation: 



{ 4|i(r), A^f,r(r) } = i2n) S,,,, 5{k k') . (6.11) 

Inserting the expansion (|2.3ip into the right-hand-side of 
Eq. (|6.10|) . we can relate the A\^\{t) mode operators to 



}(V)t 



the A'"^^ mode operators. We find: 



3.12) 



with C, = m/{ujf) + mk). Here sin^j, = l/^Jk"^ + 1 and 
cos 9 k = fc/vfi^+l, with < 0fc < TT. Density matrices 



where 



Cl^^'\)=F!,^^ir)F^^(r). 



(6.13) 
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FIG. 2: (Color online) Proper-time evolution of the matter 
components of the renormalized energy-momentum tensor for 
the one-field and two-field methods described in text. 



Particles are defined in reference to these initial states 
where a clear distinction between particles and anti- 
particles can be made. We define an average phase space 
number density Ukir) by: 



nkir) 



d2iV(r) 



dfc dr] 

and is computed using the relation: 

(+)t/ 



(6.14) 



n,(r)(2^)5(fc-fc') = (i[,;^^^(r)4S(r)) . (6.15) 



Inserting (|6.12p into (|6.15p . and using 



( A[^^ ^ ) = <5a,_ <5a'.- (2^) S{k - k') , (6.16) 



we find: 

nkir) 



\c\;-\r)?-\F^S\r)Ft\r)\^ 



(6.17) 



We see immediately that r7,fc(ro) = at t = tq. 

We note that in the one-field method the current is 
automatically zero at t = tq: Eq. (j3.7p with 



TTfc 



(6.18) 



leads to a zero current because the integrand in Eq. (|3.7p 
is odd in Tr^. Furthermore, one of the subtleties of the 
one-field method is that the zero-current point is an un- 
stable equilibrium point. This is most easily seen from 
the equation of motion, Eq. ()2.26p . of the polarization 
vector. For t — tq, we find that 



9.pW(ro) = 2K,(ro)xPi+)(r„) 

= 2Kfc(To) xKfe(ro)/cjo = 0. 

However the second derivative is not zero: 
2m2 



1 



( fc — eEo I vr? 



(6.19) 



(6.20) 



B. Two field method 

Here we start from solutions of the second-order Dirac 
equation. Writing the spinor Fk{T) in the form: 



(6.21) 



from Dirac's Eq. (|2.22p . we can find a second-order equa- 
tion for cither the upper or lower component: 



{a,2+c.2(r)-^.^fe(r)}4+;(r)==0 



(6.22) 



where s = ±1 designates the upper or lower component. 
A paramctrization of these mode functions of the form: 



X exp'^ 



A 



(+) 

k.s 



dr' 



leads to a second-order nonlinear equation for il^. ^ 
given by: 



(6.23) 



Ml 



in. 



1 STTk 

2 



(6.24) 
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two-field method ■ 




10 20 30 40 50 60 70 80 90 100 

T 

FIG. 3: (Color online) Proper-time evolution of the parti- 
cle density, dA''/d77 and proper time evolution of the ratio 
t£{t) /[AN /Arj\ for the one-field and two-field methods. 



Here, and in the following, we suppress the dependencies 
on T, fc, and the positive energy superscript. Solutions 
of the nonlinear equation (|6.24p for fig, subject to initial 
conditions given below, completely determine fs- Once 
we find fs, we can get the other Dirac component from 
Dirac's equation: 

f-s = -{idr + snk)fs = —fs, (6.25) 
771 m 

where z'l^^ (r) is given by: 

Zs = Xs+iYs = », + .s^fc - z , (6.26) 

The normalization requirement: l/sP ^ ^ means 
that: 



\f~s\' 



(6.27) 



which fixes the normalization factor As- It is an easy 
matter now to get all the terms of the density matrix ps , 
and we find: 



Pi,, = 



2Xs 


2 


+ \Zs 






2Ys 




rn? 


+ \Zs 






-\Zs 






+ \Zs 





(6.28a) 
(6.28b) 
(6.28c) 



We are now in a position to carry out an adiabatic 
expansion of the nonlinear equation (|6.24p . We again 
count derivatives with respect to r by putting: dr ^ e 9r, 
and expand 



(6.29) 



Inserting this into (|6.24p and inverting the equation gives 
0, from which we find: 



= u and ni^'^ 



2 ^ 



2uj 

From this we find that 
Zs^Xs + iYs 



5 ^2 

4 uj 



I {u) + STTk) 



(6.30) 
(6.31) 



[uj - silk) 



1 + ie 



(2) 



So from our general expressions (|6.28p , it is easy to show 
that: 



P2:s ~ em 



■2 ( ^T^l 



1 TTfc %k 5 Til ill 



4 w5 



2u^ 



(6.32) 



A. 



STTk_ 



4 UJ^ 



8 Uj'' 



For s = 1, Eqs. (|6.32p are in agreement with Eqs. (|4.4p . 
So to second adiabatic order Pi-s is independent of s, but 
P2-S and P^-s change sign with s. 

To specify the initial conditions for second-order non- 
linear Eq. (|6.24p at T = To = 1/m one needs two initial 
conditions. Since the vacuum state is not unique when 
particles are being produced, one usually chooses some 
approximate adiabatic vacuum state of given order as 
discussed in Ref. [ll|. The authors in Rcf. chose the 
first-order adiabatic conditions as 

(^") =^Q = m Vk^ + 1 , (6.33a) 

^kj(ja) =ujo = m ^ ^ , (6.33b) 

where Eq = eEQ/m?. The initial conditions are indepen- 
dent of s. 

Using Eq. (|6.32p . we obtain that for each value of s 
this choice of initial conditions at r = tq will lead to 
a non- vanishing current, Js(to). However, if we average 
over the two sets of solutions s = ± and chooses for the 
Maxwell equation: 



drE{T) 



+ 00 



[dTT,] Pi+'(7r,,r) + P3\-J(^.,r) 



(6.34) 

then the renormalizcd Maxwell equation will start with 
a zero value for the current. 
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FIG. 4: (Color online) Proper-time evolution of the momentum dependent particle density distribution, n,rfc defined in 
Eq. (|6.f4p . showing the oscillation of the centroid of the particle-density distribution between positive and negative values 
of TTfc. Here, we show results for the one-field method. Results for the two-field initial conditions scenario (not shown) are very 
similar, as it is to be expected from the results illustrated in Fig. O 



VII. NUMERICAL RESULTS 

We have performed numerical calculations for both 
sets of initial conditions described above. We employed 
a fourth-order Runge-Kutta method to solve the cou- 
pled Dirac equation and backreaction problem. The 
/c-momentum variable, which is dimensionless, was dis- 
cretized on a nonuniform piece- wise momentum grid with 
a cutoff at fc = Afe. Wc found that a value of r; 200 
was necessary to obtain numerical results insensitive with 
respect to the cutoff. For the purpose of calculating 
the subtracted values of the current J(t), matter energy 
e(r), matter pressure p{t), and fermion particle density 
dA'^(T)/d77, we needed to compute the momentum inte- 
grals with respect to the variable, tt^ rather than k. The 
corresponding momentum cutoff in TTfc-space was chosen 
to be 20% greater than TmaxAfc to allow for possible very 
large values of A{t), which is unknown at the beginning 



of the calculation. The momentum integrals in Tr^-space 
were performed using a Chebyshev integration method 
with spectral convergence [l^ . Using the procedure out- 
lined here, we found that approximately 8000 mode func- 
tions were necessary to obtain a converged numerical 
result. The conservation of the energy- momentum ten- 
sor, sec Eq. (|5.8p . served as a numerical test: we found 
that the rcnormalizcd energy-momentum tensor was con- 
served within machine precision. 

For the purpose of this comparison, we took: m = 1, 
e = 1, To = 1/m = 1, A{tq) = 0, and E{tq) = 4. These 
strong-field initial conditions have been shown to produce 
sufficient fermion pairs at t = tq for plasma oscillations 
to take place. In Fig. [U we show the proper-time evo- 
lution of the electromagnetic field, A{t), electric field, 
E{t), and current, J(t), for the one-field and two-field 
methods described in text. The components of the mat- 
ter part of the energy-momentum tensor, £(r) and p{t), 
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for the two simulations are shown in Fig. [2l Finally, the 
proper-time evolution of the particle density, dN/drj, de- 
fined in Eq. (|6.14p . is given in Fig. [3l For both methods, 
the ratio Te{T)/[dN/dri] is seen to oscillate around the 
numerical value of 1, consistent with the hydrodynami- 
cal picture, as explained in Ref.0. We notice that the two 
sets of solutions arc almost identical at short and inter- 
mediate times. The two solutions become out of phase at 
late times due to the slightly different initial conditions. 
However, in the real problem we expect that interactions 
between the fermions would eliminate these oscillations. 

The proper-time evolution of the momentum- 
dependent particle-density distribution, n^^. , correspond- 
ing to the choice of initial conditions in the one-field 
method, is shown in Fig. [H We note that the centroid of 
the particle-density distribution oscillates between pos- 
itive and negative values of tt/j. The oscillation of the 
number density is a result of the current oscillating in 
sign, the current in momentum space being related to 
the number density times the velocity of light. This ef- 
fect is also seen classically when two infinite oppositely 
charged parallel plates initially a finite distance apart are 
released and allowed to pass through one another. In that 
case both the current and electric field oscillate in an an- 
alytically derivable manner [l^. Results for the case of 
the two-field method (not shown) are very similar, as it 
is to be expected from the results depicted in Fig. [3] (see 
also Ref.[li). 



conditions which are consistent with having the fermions 
in a "vacuum state." We conclude that the two meth- 
ods of starting out the calculation produce essentially the 
same answer. Based on our numerical simulations, there 
seems to be little reason theoretically or otherwise to use 
the two-field method discussed previously in Ref. 1^, as 
it doubles the storage requirements and computational 
time. This is important for our forthcoming studies of 
fermion particle production with backreaction in QED 
and QCD. 

We emphasize here that in the case of the squared 
Dirac equation (two-field method) there are two indepen- 
dent solutions of the second-order differential equation 
for the mode functions, each of which provides a basis 
for two different fermi fields. In order to make compat- 
ible the physical requirement that the initial current is 
zero with the initial choice that the fermions were ini- 
tially chosen to be a first-order adiabatic vacuum state, 
the authors of Ref. simply averaged these two solutions 
to produce a current which was zero at r = tq. So dou- 
bling the number of fermi fields allows one to produces 
consistent initial conditions if we define the current by 
averaging over the two sets of solutions. By staying with 
the original first-order Dirac equation, in the one-field 
method we were able to satisfy the initial condition of 
zero current by choosing a slightly cruder initial state for 
the fermion fields. This choice, however, reduces by half 
the size and duration of the calculation. 



VIII. CONCLUSIONS 

To conclude, in this paper we report an initial- 
conditions sensitivity study for the problem of pair pro- 
duction of fermions coupled to a "classical" electromag- 
netic field with backreaction in (1+1) boost-invariant co- 
ordinates. We discuss two methods of choosing the initial 
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